Shock wave theory for rupture of rubber 
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This article presents a theory for the rupture of rubber. Unlike conventional cracks, ruptures in rubber travel 
faster than the speed of sound, and consist in two oblique shocks that meet at a point. Physical features of rubber 
needed for this phenomenon include Kelvin dissipation and an increase of toughness as rubber retracts. There 
are three levels of theoretical description: an approximate continuum theory, an exact analytical solution of a 
slightly simplified discrete problem, and numerical solution of realistic and fully nonlinear equations of motion. 

PACS numbers: 62.20,62.30.+d,43.25.Cbd 
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Introduction- Rapidly moving cracks in brittle materials 
under tension have a number of common characteristics. They 
cannot move faster than the shear wave speedflllll, and of- 
ten exhibit a limiting velocity around half that value because 
of instabilities of the crack tip|4]. Stresses near the tip rise 
in a universal singularity as 1/ 'y/r. In this Letter, I show that 
ruptures in rubber are different. They are supersonic. There 
is stress enhancement but no stress singularity near their tips. 
They constitute a new sort of failure mode that combines char- 
acteristics of shocks and cracks. 

The motivation for this study comes from experiments 
showing that cracks in rubber travel faster than the shear wave 
speed, and that the tip has a wedge-like shape resembling a 
shock|5]. Planar shock fronts in rubber were previously ob- 
served by Kolsky|6|. It has not been clear how to interpret 
the experiments because the large nonlinearities of rubber in- 
validate immediate comparison with the customary theory of 
linear elastic fracture mechanics. In particular, one assump- 
tion of conventional fracture mechanics is that material ahead 
of a crack tip is strained by a vanishingly small amount, while 
in popping rubber the strains are several hundred percent. 

Intersonic tensile cracks have been observed in numerics 
of Buehler, Abraham and Gao[7]. In their calculations, this 
behavior is produced by a rise in sound speed near the crack 
tip. Here the mechanism is different; there is no rise in sound 
speed. Instead, two other physical ingredients work together 
both in numerical simulations and in analytical calculations 
to reproduce the basic experimental observations. First and 
most important, the equation of motion for rubber includes 
dissipation of the Kelvin form; Langer|8] has observed that 
such terms may permit supersonic motion. Second, the rubber 
must be able to sustain larger stresses when it is relaxed along 
one axis than when it is stretched equally in all directions. 

Continuum Theory of Rubber- Strains in rubber are sev- 
eral hundred percent at rupture and one must use nonlinear 
elastic theory to describe the situation. Sound speeds in rubber 
are adequately described|5] by one of the most familiar free 
energies for non-linear elastic solids, the one due to Mooney 
and RivlinU [kJ Hill . For this free energy, define the La- 
grangean strain tensor| 12] 
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Here u(r) describes the distance from the origin of a mass 
point that was located at r before the rubber was stretched 
up. From this strain tensor one can define three rotation- 
ally invariant quantities, which are lf D = TrE, I^ = 



E 1p 



and Il D = dct E. The Mooney- 
Rivlin theory says that the free energy density of rubber is 

U/p = w = a(lf D +bll D ), (2) 



where U has units of energy per volume, p is mass density, a 
is a constant with units of velocity squared, and b is dimen- 
sionless. For a thin sheet of rubber, one can replace the three- 
dimensional theory by an effective two-dimensional one, us- 
ing the facts that rubber is highly incompressible|9|, and that 
one can neglect all the components of the strain tensor E az ex- 
cept for E zz . In two dimensions one has only two invariants, 



h = E x 



E„ 



E E 

^xx^yy 



El 



and using incompressibility to solve for E zz one finds 
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Thus one obtains an effective two-dimensional Mooney- 
Rivlin theory 



w{h,I 2 ) =a(h + bl 2 + E zz (l + bl 2 )) ■ 



(5) 



For large strains, E zz becomes negligibly small compared to 
E xx or E yy . However, as rubber relaxes to equilibrium, the 
terms proportional to E zz become important. They are what 
ensure that u = f is a minimum energy state. 

For studying the rupture of rubber, the energy density in 
Eq. 10 is both too simple and too complicated. It is too 
simple because it does not account for the fact that when 
rubber is stretched enough, the polymers pull apart and the 
force between adjacent regions drops irreversibly to zero. It 
is too complicated because the terms involving I2 and E zz 
produce nonlinear equations of motion that are impossible to 
solve analytically. Therefore, to analyze the problem, I will 
pursue two different routes. First, I will discuss numerical 
routines that supplement Eq. with information about rup- 
ture, toughening, and dissipation, and produce supersonic so- 
lutions. Second, I will isolate from Eq. (JSJi terms that are suf- 
ficient to produce good agreement with numerics and exper- 
iment, while simplifying matters enough to permit analytical 
solution. 
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Numerical System- To study rubber rupture numerically, 
consider a collection of mass points Uj whose equilibrium lo- 
cations lie on a triangular lattice, and that are connected with 
bonds to nearest neighbors. Take the lattice spacing of the un- 
stretched configuration to be A. For numerical representation 
of the strain invariants, let ulj = Uj — Hi, let n(i) refer to the 
nearest neighbors of i, and define 




if Uij < Xj 
else 



(6a) 




Figure 1: Solution of Eq. with experimental values of a — 501 
< A/ m 2 /s 2 , b — .106, Kelvin dissipation /J — 3, and rupture extension 

Gi = — y ^/, 9 a9 \ 2 (6b) A/= 5.5. Before rupture, the system is stretched vertically by a fac- 

tor of X y — 3.2, and horizontally by a factor of X x = 2.1. The 
system has been allowed to run for 500 time units, by which time 
it is approaching a steady state, apart from bunching up of material 
1 ^— \ , . . . /_ _ 2 \ 2 le - as at collides with rigid supports at top and bottom of the system. 

* = 27 h(Uij)ti{Uik) (Uij ■• u ik + 2 A J , (6c) The rupture is able to run as long as needed to the right by pasting 

j^Lk£n{%) new material on the right and discarding it from the left. Note the 

wedge-like shape of the rupture. The (Lagrangean) speed of shear 
waves ahead of the rupture is c = 21.8 m/s, and the rupture travels 
and h(u) = 1/(1 + e («-«=)/«*). (gd) at a speed v = 24.9 m/s. The system is 200 rows high with 70,000 

particles, and in its unstretched configuration is twice as wide as it is 
From these numerical quantities, one can form representa- tall, 
tions of the strain invariants as follows: 



I\ = ^/A 2 (7a) 



4 = (9/8) (G< - Hi + 4) /A 4 . (7b) 
and finally construct the energy from 

U = Y,mw(Il,P 2 ), (8) 

t 

where m is the mass in a unit cell, and the energy density w is 
given by Eq. (|5j- The quantities in (|6} are chosen according 
to two ideas. First, they are designed so that when all bonds 
at a node are shorter than a critical failure extension Xf, in the 
continuum approximation Eqs. Q reproduce the strain invari- 
ants in Eq. l|3}- Second, they are designed so that when bonds 
are stretched to an extension greater than A t , they break. For 
the three-body term in Eq. drjcl . it is necessary to introduce 
a soft cutoff through the function h described in Eq. J6di . in 
which u s is a parameter on the order of 0. 1 that sets the scale 
over which contributions to the three-body term drop to zero. 

Figure^shows an image of a steady state obtained by solv- 
ing dynamical equations that follow from Eq. The precise 
equation of motion includes dissipation of the Kelvin form, 
and is 

muf = -dU/du* + Y, ~ ««)■ (9) 

jen(i) 

One final rule is employed. Whenever some bond drops 
to a length less than 1.5 A, the failure extension A/ for the 



remaining bonds attached to nodes i and j increases. Without 
some rule of this type, the back faces of the crack disintegrate. 
Essentially, the back faces of the rupture act like a string under 
tension pulling bonds at the tip apart, and they must be able to 
sustain tensions sufficient to do so; for details, see Eq. 1161 

Numerical solutions of Eq. agree acceptably with ex- 
periment. I have tried to determine which terms in it are re- 
ally needed. Progressively stripping elements from l|9} I found 
what is most likely the simplest set of equations supporting su- 
personic solutions. These explain the nature of the solutions, 
and the conditions under which they arise. 

Neo-Hookean Continuum Theory- Experimentally! 5|, 
the dimensionless parameter b in Eq. Q is .106, so in a first 
theoretical account one can set b — 0. For strains large enough 
also to neglect E zz , Eq. reduces (up to an additive con- 
stant) to the Neo-Hookean energy density 




(10) 



The equation of motion that follows from this energy is (for 

a = x or y ) 

u a = c 2 V 2 u a . (11) 

Despite the fact that the equation of motion Eq. M \\ is an or- 
dinary wave equation, it describes large extensions. The most 
undesirable feature of this theory is that its ground state con- 
sists in material that has collapsed down to a point; this results 
from dropping the terms proportional to E zz from Eq. (|5jl. 
However, as shown in Figure|2]rupture speeds are essentially 
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unaffected by the presence of these terms. The theory has 
the great advantage that it can be solved exactly. For crack- 
like solutions, with Lagrangean variables (x, y) = ?"one has 
for y — and x < the boundary condition du y /dy = 0. 
Neither this boundary condition nor the equations of motion 
couple u x and u y ; therefore, the equations support solutions 
where u x — X x x does not change in time, and the motion of 
the mass points is purely vertical. These solutions are iden- 
tical to the solutions for a crack in anti-plane shear ! 13ll . as 
recorded for example in Ref. |1], p. 356. The static solu- 
tion has a parabolic tip. Steady states moving at velocity v, 
are identical to the static solution, but are Lorenz contracted 
in the direction of motion by a factor of yl — v 2 /c 2 . As the 
crack speed v approaches the wave speed c, the tip becomes 
increasingly blunt. 

The wave speed c has been thought the upper speed limit 
for crack-like solutions of Eq. il It . However, supersonic so- 
lutions are possible if one adds Kelvin dissipation correspond- 
ing to the rightmost term in Eq. (|9jl to obtain for a steady state 
moving at velocity v, 



2 d 2 u 2r7 2 
V = C V U 

ox z 



vc 2 /3V 



, du 
dx 



(12a) 



The variable u in this equation is the vertical motion of mass 
points u y ; the horizontal locations of all mass points remain 
fixed at u x = X x x, so there is no need to keep track of them 
further. Supplement Eq. (I12a> with the boundary conditions 



du d 2 u 
dy dxdy 



for x < 0; u = 



y 



oo. 



for x > 0. 

(12b) 

(12c) 
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/ Continuum Approximation 




□ Direct integration, N = 200 




— Discrete Solution, N = 200 




Discrete Solution, N = 2000 
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Figure 2: Comparison of theory, experiment, and numerics for 
rubber rupture velocities. Experimental velocities are scaled by 



with c —22 m/s, while the vertical extension X y is scaled by 




•, . — Ai)/3. The continuum approximation is given in Eq. 

I15bl . Direct integration of Eq. |9j is carried out in triangular lat- 
tices N =200 rows high in the Neo-Hookean limit where b — 0, 
with Kelvin dissipation /3 = 3, and retaining E zz as in Eq. 0. The 
Discrete Analytical Solution is an exact solution of the same system 
using the Wiener-Hopf technique, with the three differences. First, 
E zz is neglected in the analytical solution. Second, the analytical 
system is infinitely long in the horizontal direction, while the numer- 
ical system is finite. Third, in the numerical system there is a brief 
time when only one of two crack-line bonds has snapped, and hori- 
zontal forces on crack-line atoms do not balance to zero, while in the 
analytical solutions, all forces in the horizontal direction are ignored. 
Analytical solutions for systems both 200 and 2000 rows high are 
displayed to show how the continuum limit is achieved. Experimen- 
tal results courtesy of Paul Petersan and Robert Deegan. 



Solutions of Eq. El can be obtained with the Wiener-Hopf 
technique 1 14]. One has the following results for the upper 
face of the rupture where y = + and x<0: 



dx 



dx' 



\ y e 



x/v(3 
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Therefore, the slope a of the back face of the rupture seen in 
the lab is 



-A,, 



1 



(15a) 



This is the slope of a shock cone trailing an object traveling 
at speed v > c in a medium of wave speed c. Note that the 
velocities v and c are measured in a Lagrangean reference 
frame described by variables x and y. Horizontal and verti- 
cal speeds measured in the laboratory are larger by factors of 



X x and X y respectively; this geometrical fact accounts for the 
factor X y /X x in dl5a> . The vertical strain at the origin is ob- 
tained by setting x = in Eq. (I14> . One obtains a simple but 
approximate prediction for rupture speed by checking when 
bonds angled at 60° in a triangular lattice reach their breaking 
point Xf : 
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(15b) 
(15c) 



In order to compare with experiment, there is a single free 
parameter to fix, which is the breaking point A / . Figure [2] 
shows a comparison of the predictions from Eqs. ^]with ex- 
perimental and numerical data, using Xf — 5.5. 

An additional interesting quantity to check is the distance 
squared between horizontal mass points behind the rupture. It 
is 



K + 



v 2 /c 2 



-A f + -A x 



X , 



(16) 
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This quantity exceeds Xj for characteristic values of A/, X x , 
and X y , which explains why it is necessary for A/ to increase 
behind the rupture if the back surface is not to disintegrate. 

Neo-Hookean Discrete Theory- Not only can the contin- 
uum Neo-Hookean theory be solved, but the discrete theory, 
Eq. can also be solved exactly, provided in Eq. (0 one 
sets 6 = and E zz = 0. The solution involves the application 
of methods described in Refs. 1151 llo. Il7ll . and details will 
be presented elsewhere. Figure [2] shows exact solutions for 
rupture speeds in systems 200 and 2000 rows high compared 
both with direct integration of the equations of motion and 
experiment. In addition to removing discrepancies between 
the very simple results in Eqs. J15i and numerics, solving the 
discrete model explains the conditions under which one gets 
supersonic or subsonic solutions for cracks in tension. 

The basic result is this: including dissipation through (3 in 
the equation of motion introduces a length scale (3c into the 
problem. The behavior of cracks hinges on the ratio of (3c 
to the lattice spacing A. When (3cj A is much less than one, 
cracks behave as in conventional fracture mechanics, and their 
speed is limited from above by c, except within a very narrow 
window of strains where all bonds in the system ahead of the 
crack approach their breaking point. As (3c/ A approaches and 
exceeds one, dissipation progressively destroys the stress sin- 
gularity around conventional crack solutions, but at the same 
time it permits the appearance of supersonic solutions. Note in 
Eq. dl5b> that rupture speed is determined by vertical exten- 
sion X y , rather than by the total energy stored ahead of the 



crack tip as in conventional fracture mechanics. Exact so- 
lution of the discrete Neo-Hookean theory shows that dl5bt 
is not completely accurate, but its scaling properties are cor- 
rect. One sees in Figure [2] that the relation between rupture 
velocity and system extension X y has essentially reached the 
macroscopic limit for systems 200 rows high and velocities v 
above 1.05c. The macroscopic limit is subtle near v = c, since 
solutions with speeds above and below c scale differently as 
system size goes to infinity. 

Establishing the existence of supersonic ruptures in tension 
opens up many possibilities for future work. The supersonic 
ruptures in experiment begin to oscillate once exceeds a 
critical value. The numerical and analytical tools provided 
here should provide an appropriate starting point for studying 
the oscillations. Finally, it would be interesting to know if 
there are materials different from rubber that meet the condi- 
tions needed to sustain supersonic ruptures. 



Acknowledgments 

I am indebted to Jim Rice for pointing out, in a lengthy 
email, that it would be profitable to study this problem with 
the Neo-Hookean theory. I have had many discussions about 
the physics with Robert Deegan, Paul Petersan, and Harry 
Swinney. Financial support from the National Science Foun- 
dation through DMR-0401766 is gratefully acknowledged. 



[1] L. B. Freund, Dynamic Fracture Mechanics (Cambridge Uni- 
versity Press, Cambridge, 1990). 

[2] K. B. Broberg, Cracks and Fracture (Academic Press, San 
Diego, 1999). 

[3] K. Ravi-Chandar, in Comprehensive Structural Integrity, edited 
by I. Milne, R. O. Ritchie, and B. Karihaloo (Elsevier, 2003), 
vol. 2, chap. 5. 

[4] J. Fineberg and M. Marder, Physics Reports 313, 1 (1999). 
[5] P. J. Petersan, R. D. Deegan, M. Marder, and H. L. Swinney, 

PRL 015504 (2004). 
[6] H. Kolsky, Nature 224, 1301 (1969). 

[7] M. J. Buehler, F. F. Abraham, and H. Gao, Nature 426, 141 
(2003). 

[8] J. S. Langer, Physical Review A 46, 3123 (1992). 
[9] L. R. G. Treloar, The Physics of Rubber Elasticity (Oxford Uni- 
versity Press, Oxford, 1975), 3rd ed. 



[10] M. Mooney, Journal of Applied Physics 11, 582 (1940). 

[11] R. Rivlin, Philosophical Transactions A 241, 379 (1948). 

[12] A. C. Eringen and E. S. Suhubi, Elastodynamics, vol. 1 (Aca- 
demic Press, New York, 1974). 

[13] W. W. Klingbeil and R. T. Shield, Zeitschrift fur angewandte 
Mathematik und Physik 17, 281 (1966). 

[14] B. Noble, Methods Based on the Wiener-Hopf Technique for 
the Solution of Partial Differential Equations (Pergamon, New 
York, 1958). 

[15] L. I. Slepyan, Models and Phenomena in Fracture Mechanics 

(Springer, Berlin, 2002). 
[16] S. I. Heizler, D. A. Kessler, and H. Levine, PRE 66, 016126/1 

(2002). 

[17] M. Marder and S. Gross, Journal of the Mechanics and Physics 
of Solids 43, 1 (1995). 



